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Abstract: We investigate neutrino oscillations from core-collapse supernovae that pro- 
duce magnetohydrodynamic (MHD) explosions. By calculating numerically the flavor con- 
version of neutrinos in the highly non-spherical envelope, we study how the explosion 
anisotropy has impacts on the emergent neutrino spectra through the Mikheyev-Smirnov- 
Wolfenstein effect. In the case of the inverted mass hierarchy with a relatively large #13 
(sin 2 2#i3 > 1CP 3 ), we show that survival probabilities of v e and v e seen from the rotational 
axis of the MHD supernovae (i.e., polar direction), can be significantly different from those 
along the equatorial direction. The event numbers of v e observed from the polar direction 
are predicted to show steepest decrease, reflecting the passage of the magneto-driven shock 
to the so-called high-resonance regions. Furthermore we point out that such a shock effect, 
depending on the original neutrino spectra, appears also for the low-resonance regions, 
which could lead to a noticeable decrease in the v e signals. This reflects a unique nature 
of the magnetic explosion featuring a very early shock-arrival to the resonance regions, 
which is in sharp contrast to the neutrino-driven delayed supernova models. Our results 
suggest that the two features in the v e and v e signals, if visible to the Super-Kamiokande 
for a Galactic supernova, could mark an observational signature of the magnetically driven 
explosions, presumably linked to the formation of magnetars and/or long-duration gamma- 
ray bursts. 



Keywords: 



neutrino conversion;matter effects;supernova;magnetic fields;rotation. 



Contents 



^ Introduction |l| 

[2L Numerical method |3| 

|2.1| MHD driven supernova explosion model || 

Neutrino oscillation calculations 



||. Result @ 
|3.l| Survival probability of v e @ 
3.2 Neutrino spectra of v e I I 



|3.3| Expected event number of v e at the Super-Kamiokande detector 11 
Expected event number of v e at the Super-Kamiokande detector 13 



Discussion 15 
[4.1| Dependence of #13 on neutrinos 15 
Difference of the original neutrino spectrum 16 



IHl Summary 19 



1. Introduction 

Current estimates of Galactic core-collapse supernovae rates predict one core-collapse su- 
pernova event every ~ 40±10 year (e.g., Q). When a massive star (> 1OM ) |^] undergoes 
a core-collapse supernova explosion in our Galactic center, copious numbers of neutrinos 
are produced, some of which may be detected on the earth. Such supernova neutrinos will 
carry valuable information from deep inside the core. In fact, the detection of neutrinos 
from SN1987A (albeit in the Large Magellanic Cloud) paved the way for Neutrino As- 
tronomy, an alternative to conventional astronomy by electromagnetic waves [||, ||]. Even 
though neutrino events from SN1987A were just two dozens, they have been studied exten- 
sively and allowed us to have a confidence that the basic picture of core-collapse supernova 
is correct (|||, || see 0] for a recent review). Over the last decades, significant progress has 
been made in a ground-based large water Cherenkov detector as Super-Kamiokande (SK) 
Q and also in a liquid scintillator detector as KamL AND || . If a supernova occurs in our 



Galactic center (~ 10 kpc), about 10,000 v e events are estimated to be detected by SK [10 



(see also [11]). Those successful neutrino detections are important not only to study the 
supernova physics but also to unveil the nature of neutrinos itself such as the neutrino os- 
cillation parameters the mass hierarchy, and neutrino mass [12, 12, 14, 15, [If], 17, [l~8||(e.g., 
[IS] for a recent review). 
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The neutrino flux we receive on the Earth is sensitive to the profile of the matter 
that the neutrinos encounter on their path. The supernova neutrinos interact with elec- 
trons when the neutrinos propagate through stellar matter via the well-known Mikheyev- 
bmirnov- Wolfenstein (MSW) effect pL M, EH1 p|. Such effects have been extensively 



investigated so far from various points of view, with a focus such as on the progenitor 
dependence of the early neutrino burst [23, 24] and on the Earth matter effects (e.g., 

mm)- 

The neutrino conversion efficiency via the MSW effect depends sharply on the den- 
sity/electron fraction gradients, thus sensitive to the discontinuity produced by the passage 
of the supernova shock. It is rather recently that such shock effects have been focused on 
p8|| (sec |HJ for a recent review). The time dependence of the events showing a 



BE 

decrease of the average energy of v e in the case of normal mass hierarchy (or v e in the 
case of inverted mass hierarchy), monitors the evolution of the density profile like a to- 
mography, thus could provide a powerful test of the mixing angle and the mass hierarchy 
(e.g., [17, [K], |3l|]). Very recently the flavor conversions driven by neutrino self- interactions 
are attracting great attention, because they can induce dramatic observable effects such as 
by the spectral split or swap (e.g., 33, 34, 35] and see references therein). They are 
predicted to emerge as a distinct feature in their energy spectra (see |2^, for reviews 
of the rapidly growing research field and collective references therein). The resonant spin- 

p9| and references 



flavor conversion has been also studied both analytically (e.g., [37, [ 
therein) and numerically (e.g., [40, 41 1), which induces flavor conversions between neutri- 
nos and antineutrinos in magnetized supernova envelopes. Those important ingredients 
related to the flavor conversions in the supernova environment have been studied often one 
by one in each study without putting all the effects together, possibly in order to highlight 
the new ingredient. In this sense, all the studies mentioned above should be regarded as 
complimentary towards the precise predictions of supernova neutrinos. 

Here it should be mentioned that most of those rich phenomenology of supernova neu- 
trinos have been based on the spherically symmetric models of core-collapse supernovae 



p8| . On the other hand, there are accumulating observations indicating that 
core-collapse supernovae are globally aspherical commonly (e.g., |4^]). Pushed by them, 
various mechanisms have been explored thus far by supernova modelers to understand 
the central engine, such as the neutrino-heating mechanism aided by convection and the 
standing accretion shock instability (e.g., [43, 44 1 and references therein), magnetohydro- 
dynamic (MHD) mechanism relying on the extraction of the rotational energy of rapidly 
rotating protoneutron stars via magnetic fields (e.g., [45, 4£, 47], see for a recent re- 
view and collective references therein) , and the acoustic mechanism relying on the acoustic 
energy deposition via oscillating protoneutron stars [^]. Apart from simple parametriza- 
tion to mimic anisotropy and stochasticity of the shock (e.g., [49, 50|), there have been a 
very few studies focusing how those anisotropies obtained by the recent supernova simu- 
lations could affect the neutrino oscillations ([30, 51]). This is mainly due to the lack of 
multidimensional supernova models, which are generally too computationally expensive to 
continue the simulations till the shock waves propagate outward until they affect neutrino 
transformations. 
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Figure 1: Snapshots showing MHD explosions of core-collapse supernovae at 0.8, 2.0 and 4.3 sec 
after core bounce from left to right. In each figure, contour of the logarithmic density [g/cm 3 ] (left) 
and entropy per baryon [ks] (right) are shown. The unit of the horizontal and the vertical axis is 
in [cm]. Note that the difference of the length scale for each panel and that the outermost radius 
of the progenitor is ~ 3 x 10 10 [cm]. 



Here we study the neutrino oscillations in the case of the MHD explosions of core- 
collapse supernovae. In this case, the highly collimated shock pushed by the strong mag- 
netic pressure can blow up the massive stars along the rotational axis [4(], [47|. Such 



explosions are considered to precede the formations of magnetars (e.g., [45, 47, p2|), pre- 



sumably linked to the so-called collapsars (e.g., |53|, |54J, |55j), which attract great attention 



recently as a central engine of long-duration gamma-ray bursts (GRBs) (e.g., |56|]). By 



performing special relativistic simulations [47], it is recently made possible to follow the 



MHD explosions in a consistent manner, starting from the onset of gravitational collapse, 
through core-bounce, the magnetic shock-revival, till the shock-propagation to the stellar 
surface. Based on such models, we calculate numerically the flavor conversion in the highly 
non-spherical envelope through the MSW effect. For simplicity, we neglect the effects 
of neutrino self-interactions and the resonant spin-flavor conversion in this study. Even 
though far from comprehensive in this respect, we present the first discussion how the 
magneto-driven explosion anisotropy has impacts on the emergent neutrino spectra and 
the resulting event number observed by the SK for a future Galactic supernova. 

The paper opens with generalities on magneto-driven explosion models and details on 
our numerical model (section 2). In section 3, main results are shown, in which we pay 
attention to the shock-effects on the so-called high-resonance in section 3.1-3.3 and on the 
low-resonance in section 3.4, respectively. In section 4, we give a discussion focusing on 
the dependence of #13 and also of the original neutrino spectra on the obtained results. 
Summary is given in section 5 with implications of our findings. 

2. Numerical method 

2.1 MHD driven supernova explosion model 

We take time-dependent density profiles from our 2.5 dimensional MHD simulations of 



core-collapse supernovae [47], in which a 25 M presupernova model by Heger et al. [57] 



was adopted. This progenitor lacks the hydrogen and helium layer during stellar evolution, 
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which is reconciled with observations that the progenitors associated with long-duration 



gamma-ray bursts are type Ib/c supernovae (e.g., [56, p8fl). 

For our MHD model, the strong precollapse magnetic fields of ~ 10 12 G are imposed 
with the high angular velocity of ~ 72 rad/s with a quadratic cutoff at 100km in radius, 
in which the rotational energy is equal to ~1% of the gravitational energy of the iron core 
(model B12TW1.0 in |17|]). Such a combination of rapid rotation and strong magnetic 
field, although strong differential rotation is assumed here, have been often employed in 



the context of MHD supernova models jig, 47, 6C].In those models, the explosion proceeds 



by the strong magnetic pressure amplified by the differential rotation near the surface of 
the protoneutron star. 

One prominent feature of the MHD models is a high degree of the explosion asphericity. 
Figure |l| shows several snapshots featuring typical hydrodynamics of the model, from near 
core-bounce (0.8 s, left), during the shock-propagation (2.0 s, middle), till near the shock 
break-out from the star (4.3 s, right), in which time is measured from the epoch of core 
bounce. It can be seen that the strong shock propagates outwards with time along the 
rotational axis. On the other hand, the density profile hardly changes in the equatorial 
direction, which is a generic feature of the MHD explosions observed in recent numerical 



simulations p5j |47|, |60| 



As well known, the flavor conversion through the pure-matter MSW effect occurs in 
the resonance layer, where the density is 



. , Am 2 \ /10MeV\ /0.5\ 

*- ~ x 10 g/cm {wvf) {—) (n J cos2e (21) 

where Am 2 is the mass squared difference, E u is the neutrino energy, Y e is the number 
of electrons per baryon, and 9 is the mixing angle. Since the inner supernova core is too 
dense to allow MSW resonance conversion, we focus on two resonance points in the outer 
supernova envelope. One that occurs at higher density is called the H-resonance, and the 
other, which occurs at lower density, is called the L-resonance. Am? and 9 correspond to 
Am 2 3 and #13 at the H-resonance and to Am? l2 an d #12 at the L-resonance. 

Figures |2| and ||| are evolutions of the density profiles in the polar direction every 0.4 s 
as a function of radius. Above and below horizontal lines show approximately density of 
the resonance for different neutrino energies which are 5 and 60 MeV, respectively. Blue 
lines (left) show the range of the density of the H-resonance, and sky-blue lines (right) 
show that of the L-resonance. Along the polar axis, the shock wave reaches to the H- 
resonance, ~ O(10 3 )g/cm 3 at ~ 0.5 s, and the L-resonance, ~ 0(l)g/cm 3 at ~ 1.2 s. It 
should be noted that those timescales are very early in comparison with the ones predicted 
in the neutrino-driven explosion models, typically ~ 5 s and ~ 15 s for the H- and L- 



resonances, respectively (e.g., [17, 30]). This arises from the fact that the MHD explosion 
is triggered promptly after core bounce without the shock-stall, which is in sharp contrast 
to the neutrino-driven delayed explosion models { WA |30tl). The progenitor of the MHD 
models, possibly linked to long-duration gamma-ray bursts, is more compact due to the 
mass loss of hydrogen and helium envelopes during stellar evolution , which is also the 
reason for the early shock-arrival to the resonance regions. 
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Figure 2: The density profile in the po- 
lar direction every 0.4 s as a function of ra- 
dius, (a), (b) and (c) correspond to 0.4, 
0.8 and 2.8 s, respectively. The horizontal 
blue lines show approximately density of the 
H-resonance for different neutrino energies 
which are 5 (above) and 60 MeV (below), 
respectively. The green line is an original 
density profile for 2.0 s. 
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Figure 3: Same as Figure 2, but for 1.2 
(d), 2.0(e), and 3.2 (f) s, respectively. 



Figure |] shows the density profiles in the equatorial direction every 0.4 s. The profiles 
are almost unchanged, because the MHD shock does not reach to the outer layer for the 
equatorial direction. It is noted in Figures 2 and 3 that to maximize the shock effect, a 
sharpness of the shock along the polar direction is modified from the one in the MHD 
simulations, which is made to be inevitably blunt due to the employed numerical scheme 
using an artificial viscosity to capture shocks |47|]. Green line in Figure [2] is an original 
density profile for 2.0 s, whose shock front is replaced with the sharp vertical surface passing 
through the midpoint of the blunt shock. Although our simulations are one of the state-of- 
the-art MHD simulations at present, which can encompass the long-term hydrodynamics 
starting from gravitational collapse till the shock-breakout [[l7j, this is really a limitation 
of the current MHD simulations. 

Snapshots labeled as (a), (b), (c), (d), (e) and (f) in Figures [2] and || correspond to 
0.4, 0.8 and 2.8, 1.2, 2.0 and 3.2 s after bounce, respectively, which is chosen for later 
convenience to compute the time evolutions of the flavor conversions. 

2.2 Neutrino oscillation calculations 



Along the time-dependent density profiles as shown in Figures |2j-[|, we solve numerically 
the time evolution equation for the neutrino wave functions as follows, 



d 
i— 
dt 
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(2.2) 



where AEij = Am^/2£' l/ , and Am^ are mass squared differences, E v is a neutrino energy, 
Gp is the Fermi constant, and n e is an electron number density. In the case of anti-neutrino, 
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Figure 4: The density profile in equatorial 
direction every 0.4 s as a function of radius. 
The highly collimated shock does not propa- 
gate in the equatorial direction. 
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Figure 5: The original spectra at 0.4 s after 
core bounce taken from the Lawrence Liver- 
more simulation (e.g., |l(J). The solid line 
(red line), the dashed line (green) and the 
dotted line (blue) are v e , v e and v x spectra, 
respectively. 



the sign of \/2G Fn e changes. U is the Cabibbo-Kobayashi-Masukawa (CKM) matrix, 



U 



C12C13 



" Cl2S23Sl 3 e' u 
Cl2C23Sl3e l<5 



C12C23 " 
-C12S23 



S12C13 
- Sl2S23«13e' 



•Sl2C2 3 Sl 3 e 




(2.3) 



where Sij 



-S12C23 
\ S12S23 - 

sinf^jj and Cjj = cos 9^ (i ^ j = 1,2,3) are the mixing angle. We set the CP 
violating phase, 5, equal to zero in the CKM matrix for simplicity. Solving numerically the 
above equation, we obtain neutrino survival probabilities P (P) that neutrinos emitted as 
v e [pe) at the neutrino sphere remain v e (p e ) at the surface of the star |61]. The survival 
probabilities are determined by an adiabaticity parameter 7, 

Am 2 sin 2 29 



7 



2E U cos 29 (l/n e ) (dn e /dr) 



(2.4) 



Am 2 and 9 correspond to Am 2 3 and # 13 at H-resonance, and to Amf 2 and 6q 2 at L- 
resonance, respectively. When 7 > 1, the resonance is called "adiabatic" and the conver- 
sion between mass eigenstates does not occur. On the other hand, when 7 <C 1, the reso- 
nance is called "non-adiabatic" and the mass eigenstates are completely exchanged. If the 
mass hierarchy is normal, the H- and L-resonances occur only in the neutrino sector. On the 
other hand, if the hierarchy is inverted, the H-resonance occurs in the anti-neutrino sector, 
and the L-resonance occurs in the neutrino sector. The neutrino oscillation parameters are 
taken as sin 2 26» 12 =0.84, sin 2 26» 23 = 1.00, Amf 2 =8.1x 10 _5 eV 2 and | Amf 3 |=2.2x 10~ 3 eV 2 
(e.g., the summary in [62] and references therein). For the unknown properties, we assume 
inverted mass hierarchy and sin 2 2#i3=1.0x 10~ 3 as a fiducial value in our computations. 
This is because observable features in z/ e , most accessible channel to the SK, would imply 
that the neutrino mass hierarchy is inverted and that # 13 is large. The dependence of this 



parameter on the results will be discussed in Section 4.1 



The neutrino energy spectra at the surface of the star, (j)^(E v ), are calculated by 
multiplying the survival probabilities by the original supernova neutrino spectra, (j)®(E u ) 
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(2.5) 



where 0^ = 1/4(0^ + <^, r + 0p M + ^p r ). To continue the simulations till the MHD shocks 
propagate outward until they affect neutrino transformations, the protoneutron stars are 
excised when the shocks comes out of the iron core. By doing so, the severe Courant- 
Friedrich-Levy condition in the stellar center can be relaxed, which is often taken in the 
long-term supernova simulation (e.g., [44, [|l]] and references therein). A major drawback 
by the excision is that the information of the neutrino spectra at the neutrino sphere is lost. 
Assuming that the qualitative behavior in the neutrino luminosity, such as the peaking near 
neutronization and the subsequent decay with time, is similar between the neutrino-heating 
and the MHD mechanism, we employ the results of a full-scale numerical simulation by the 
Lawrence Livermore group (e.g., as in pi]]), which predict the neutrino temperatures of 
u e , v e and v x to be 2.8, 4.0, and 7.0 MeV, respectively. The neutrino spectra at the neutrino 
sphere are approximated by Fermi-Dirac distributions and the neutrino luminosity is taken 
to decay exponentially with a timescale of ti, = 3 s (e.g., [|63]]). The total neutrino energy 
is assumed to be fixed at 2.9 x 10 53 erg that corresponds to the model of the Lawrence 
Livermore group. Figure || shows the original spectra at 0.4 s after core bounce. The 
variation of the original neutrino spectra and its effects on our results will be discussed in 
Section |J. 

We calculate expected event numbers of the supernova neutrinos at a water Cherenkov 
detector. The positron (or electron) energy spectrum is evaluated as follows, 



d 2 N 



1 



d 2 N v 



dE v 



(2.6) 



dE e dt " haL c/ 4vrd 2 dE v dt " v dE e 
where N is a detection number of neutrinos, iV t ar is a target number, rj(E e ) is a efficiency 
of the detector, E e is an energy of electron or positron, d is a distance from the supernova, 
§£f^g is a neutrino spectrum from the star calculated by the procedure mentioned above, 
and a is the corresponding cross section given in |H| . We consider the SK detector which 
is filled with 32000 ton pure water. The finite energy resolution of the detector is neglected 
here. The event number is obtained by integrating over the angular distribution of the 
events. For simplification, we take the efficiency of the detector as follows: rj(E e ) = at 
E e < 7 [MeV] and r)(E e ) = 1 at E e > 7 [MeV] |6~4j |. We assume that the supernova occurs 

10 [kpc]). Albeit important 



in our Galactic center (d = 10 [kpc]). Albeit important [^5|, 26 1, we do not consider the 
Earth matter effect to highlight first the MSW effect inside the star. 



3. Result 

3.1 Survival probability of v e 

In this section, we focus on the influence of the shock wave in the H-resonance that appears 
in the v e signatures in the inverted mass hierarchy. First, we calculate survival probability 
of u e: P(E U ), in the manner described in Section |2.2j . 
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Figure 6: The survival probability of v e in the case of (a), (b) and (c). The solid lines (red lines) 
and the dotted lines (blue lines) are for the polar and equatorial direction, respectively. Inverted 
mass hierarchy and sin 2 26>i3 = are assumed. The survival probability of the polar direction 

changes from very early stage due to the faster shock-arrival to the resonance layers than the 
neutrino-driven supernova models (See text for more details). 



Figure || is the survival probability of v e as a function of neutrino energy in the case of 
(a), (b) and (c). The solid lines (red lines) and the dotted lines (blue lines) indicate that 
of the polar direction and the equatorial direction, respectively. The survival probability 
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of the equatorial direction is almost in all the panels. Here we focus on the survival 
probability of the polar direction. First, in the panel (a), the low-energy side of the survival 
probability becomes finite. Next, in the panel (b), the survival probability becomes finite 
in all the considered energy range. Finally, in the panel (c), the high-energy side of the 
polar direction only remains finite. The region of the finite survival probability shifts from 
the low-energy side to the high-energy side with time. The finite value is close to P ~ 0.7 
that is close to the case of the complete non-adiabatic state. It is noted that the survival 
probability for v e can be given as follows (e.g., fll2|| ), 



P = P H cos 2 9 12 cos 2 flia + (1 - P H ) sin 2 9 13 , (3.1) 

and with the employed value of cos 2 6*12 = 0.7, 

P = 0.7P H cos 2 0i3 + (1 - P H ) sin 2 9 13 , (3.2) 

where Pjj is the survival probability at the H-resonance. With the shock at the H-resonance, 
Ph is close to unity, which makes P close to 0.7. 

The behaviors in Figure 6 can be understood by comparing the resonance point with 
the shock position. When the shock front reaches the region of H-resonance (see the region 
enclosed by the blue lines in Figure ||), the steep decline of the density at the shock front 
changes the resonance into non-adiabatic following Eq. ( j2.4| ). As a result, the survival 
probability starts to become finite. As the shock propagate from the high density region 
to low energy region, the survival probability becomes finite from low-energy side to high- 
energy side, since the density at the resonance point, p res , is proportional to E~ x (see 
Eq. (|2 . 1]) ) . In the equatorial direction, on the other hand, the survival probability does 
not change (blue dotted lines of Figure ^). This is simply because the density near the 
resonance along the equatorial plane hardly changes as already shown in Figure 

The adiabaticity of the conversion can be evaluated by comparing the oscillation length, 
L osc , with the scale height of electron number density, 5r as follows, 

5r 

7 = 7 — > 1 <=^> 5r > L osc , (3.3) 



^osc 



2E U 
Am 2 



( 2V2G F E u n e V . 2 

- — cos 20 + sin 20 



\ Am 2 



(3.4) 



sin 20 n e 

6r = — 7^7' ( 3 - 5 ) 

cos 20 \an e /ar\ 

When the scale height is shorter than the oscillation length, the resonance is non-adiabatic. 
Three panels in Figure show the scale height and the oscillation length of neutrinos as 
a function of radius. These panels correspond to the case of (a), (b) and (c) from top to 
bottom. The dashed lines (green lines) and the solid lines (red lines) are the oscillation 
length of anti-neutrinos which energy are 5 MeV and 60 MeV, respectively. Dotted lines 
(blue lines) are the scale height. Most steepest decline of the scale height shown in each 
panel coincides with the position of the shock front. Each horizontal line indicates the 
maximum value of the oscillation length and the minimum value of the scale height. 
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Figure 7: The dotted blue lines represent the scale height (e.g., Eq.(3.3)). The dashed green lines 
and the solid red lines are the oscillation length (Eq. (3.2)) of anti-neutrinos for 5 MeV and 60 
MeV, respectively. Each horizontal lines are the maximum of oscillation length and the minimum 
values of scale height. It is noted that the label of (13) indicates that the value depends sensitively 
on and AmJ 3 . 



In Figure 0(a), the scale height is shown to be shorter than the oscillation length of 5 
MeV at the shock front at ~ O.O6i?0 with R@ being the solar radius (compare the minimum 
of the blue line with the maximum of the green line, which is shown by the horizontal green 
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line). Therefore, the H-resonance of 5 MeV becomes non-adiabatic (the survival probability 
~ 0.7). On the other hand, the scale height is larger, albeit slightly, than the oscillation 
length of 60 MeV, making the resonance adiabatic (the survival probability ~ 0). In Figure 
|7|(b), so the scale height is smaller than the oscillation length of 5 MeV and 60 MeV at the 
shock front (~ O.li?©), the H-resonance for almost all the neutrino energies become non- 
adiabatic, making the survival probability close to be 0.7. Figure ^(c) is a direct opposition 
as for the resonance conditions to the case of (a) . The oscillation length for the high-energy 
neutrinos is longer than for the low-energy neutrinos (Eq. (3.2)). This is the reason of the 
slight increase in the survival probability for the high-energy side seen in the equatorial 
direction (Figure ||). 

3.2 Neutrino spectra of v e 

Now we move on to estimate the neutrino spectra at the surface of the star using the 
calculated survival probability in the previous section, and the original neutrino spectra in 
section 2.2. 

From Eq. (|2.5|), the spectra of v e from the star ((f>p^) is expressed as, 

4: = p<t>i + (i - p)< 

= P(^ e -<) + <. (3.6) 

The original spectra of the low-energy side are <f>% e > 4>® x (see Figure ||), thus the sign of 
(cf)p e — ^ ) becomes positive. Therefore, cf>^ increases compared with the no-shock case 
because P can be finite due to the shock-passage (e.g., section 3.1). On the other hand, 
the original spectra of the high-energy side are (j)® < 0° , making the sign of ((j)® — 4>® ) 
negative, leading to the decrease in in comparison with the no-shock case. It should 
be noted that the energy satisfying (fP De = ^ is about 23.3 MeV (See Figure^). 

Figure || is the energy spectra of v e in the case of (a), (b) and (c) from top to bottom. 
The solid lines (red lines) and the dotted lines (blue lines) are the spectra of the polar 
direction and the equatorial direction, respectively. At the panel (a), an enhancement of 
the low-energy side of the polar direction is seen. This is because the survival probability 
at the low energy side is finite as shown in the panel (a) of Figure ^. The effect of the 
shock on the neutrino spectra also shifts from the low-energy to the high-energy side as 
the shock front goes outward. In fact, the energy spectrum of the polar direction becomes 
softer in the high-energy side than for the equatorial direction as shown in Figure |8] (b) . 

3.3 Expected event number of v e at the Super-Kamiokande detector 

From Eq. (2.6), we calculate the expected event number at SK for a Galactic supernova. 
Figure |9] shows the time evolution of the event number, in which the solid line (red line) 
and the dotted line (blue line) are for the polar and the equatorial direction, respectively. 
The shock effect is clearly seen. The event number of the polar direction shows a steep 
decrease, marking the shock passage to the H-resonance layer. Moreover, we can see a 
slight enhancement for the polar direction around 0.5 sec, because of the increase of the 
low-energy neutrinos by the shock propagation that is seen in the panels (a) and (b) of 
Figure ||. 
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Figure 8: v e spectra at the surface of the star in the case of (a), (b) and (c) from top to bottom. 
We assume the inverted mass hierarchy and sin 2 2#i3 = 10 -3 . The solid lines (red lines) and the 
dotted lines (blue lines) are the spectra of polar direction and equatorial direction, respectively. 



It is noted that the decrease in the events comes mainly from the decrease of the high- 
energy neutrinos rather than the increase of the low-energy neutrinos. This is because the 
cross section of v e + p — * e + + n, main reaction for detection, is proportional to the square 
of the neutrino energy (E%). The change of the event number by the shock passage is 
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Figure 9: The expected event number of v e 
in the SK as a function of the time. The solid 
line (red line) and the the dotted line (blue 
line) is for the polar and the equatorial direc- 
tion, respectively. The supernova is assumed 
to be located at the distance of 10 kpc. The 
error bars represent the 1<7 statistical errors 
only. 
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Figure 10: Same as Figure 9, but for the 
ratio of high- to low-energy events defined in 
Eq. (3.7). 



about 36% of the event number without the shock. Since the expected events are ~ 2500 
at the sudden decrease (Figure ||), it seems to be quite possible to identify such a feature 
by the SK class detectors. Such a large number imprinting the shock effect, possibly up to 
two-orders-of magnitudes larger than the ones predicted in the neutrino-driven explosions 
(e.g., Fig 11 in |3(|), is thanks to the mentioned early shock-arrival to the resonance layer, 
peculiar for magnetic explosions. 



Here we define a ratio of the high-energy to the low-energy event numbers as in [31]: 



event number of 20 < E < 60 MeV 
' ' ' ~ event number of 5 < E < 20 MeV ' ^ ' ' 

It is shown in Figure pi] that the ratio could be a helpful tool to identify the shock effect 
more clearly. 

3.4 Expected event number of u e at the Super-Kamiokande detector 

Now we proceed to discuss the shock effect on the L-resonance (Figure ||). 

In the L-resonance region, the shock influences not v e but v e in the inverted mass 
hierarchy. We calculate the survival probability of v e% P(E U ), in the same manner described 
in Section |2.2| . As will be shown below, the qualitative behavior of v e in the L-resonance 
closely resembles to that of v e in the H-resonance. Just like the case of the H-resonance 
(section |3.1| ), the survival probability of v e is about 0.3 in absence of the shock, becomes 
0.7, when the shock reaches to the L-resonance layer. Also as in the case of u e , the shock 
effect on the survival probabilities shifts from the low- to high- energy side as the shock 
propagates outward, leading to the change in the observed spectra of v e (Figure 11). 

With these survival probabilities, we calculate the expected event number at the SK 
for a Galactic source (Figure |l2|). It can be seen that the event number of the polar 
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Figure 11: The v e spectra in the case of (d), (e) and (f) from top to bottom, 
inverted mass hierarchy and sin 2 2t9 13 = 10 -3 . 



We assume the 



direction initially increases and then decreases compared with the one for the equatorial 
direction, which is due to the increase of the low-energy neutrinos and due to decrease of 
the high-energy neutrinos, respectively. The event number of v e become much fewer than 
that of v e . This is because the main reaction for detecting v e is v e + e~ — > v e + e~, and 
the cross section is about 10~ 2 times smaller than that of v e . Note that the cross section 
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Figure 12: The expected event number of 
v e in the SK as a function of time. The solid 
line (red line) and the dotted line (blue line) 
for the polar and equatorial direction, respec- 
tively. The supernova is assumed to be lo- 
cated at the distance of 10 kpc. 



o i i i i i i i I 

0.5 1 1.5 2 2.5 3 
Time [sec] 

Figure 13: Same as Figure 12, but for the 
ratio defined in Eq. (3.7). The error bars 
represent the statistical errors only. 



of f e + e — > v e + e is largest among neutrino-electron scattering reactions (e.g., |35|]) and 
that this forward scattering reaction could be an important tool for identifying the direction 



to the supernova (|5(], 67, 68, |6J|). As in the case of the P e , R(H/L) defined in Eq.(^J) 
could be a good probe, especially when the event number is smaller here (Figure |i~3| ). The 
ratio begins to decrease near around 1 s, which might be an observable signature of the 
shock entering to the L-resonance. While the statistical errors are eye-catching compared 
to those of the H-resonance due to the small event numbers here (compare Figure |||), the 
overlap of the error bars from 1.5 to 2.5s in the high/low-ratio is fortunately small, which 
could mark a marginally distinguishable feature. 

To our knowledge, the possibility of observing the shock feature in the L-resonance 
has been never claimed before. This is because in the conventional delayed explosion 
models of core-collapse supernovae, the neutrino oscillation in L-resonance is thought to 
be minor because the original neutrino luminosity decreases to be so small until the shock 
wave reaches to the L-resonance. Due to the small cross section of f e , the detection of 
the shock effect should not be easy even for the case of magnetic explosions considered 
here. However we speculate the detection to be more promising when the proposed next- 
generation megaton-class detectors such as the Hyper-Kamiokande (~0.5-Mton) |7(| and 
the Deep-TITAND (~5-Mton) [[n|] are on-line, possibly making the event number up to 
two-orders-of-magnitudes higher than for the SK. 



4. Discussion 



4.1 Dependence of #13 on neutrinos 

First of all, we discuss the dependence of #13 on the expected event number. Considering 
the uncertainty of #13, we calculate the event number changing sin 2 2^13 = 10~ 2 and 10 -5 
in addition to the fiducial value of 10 -3 . 
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Figure 14: The expected event number of v e are shown. The expected event number in the 
left panel are calculated with sin 2 2#i3 = 1CP 2 and that in the right panel are calculated with 
sin 2 2#i3 = 10~ 5 . The solid lines (red lines) and the dotted lines (blue lines) are the event number 
in the polar direction and the equatorial direction, respectively. 



The two panels of Figure |TJ show the event number of v e for sin 2 2#i3 = 1CP 2 (left) 
and 10 -5 (right). It can be seen that the event number of sin 2 2#i3 = 10 -2 is larger than 
that of sin 2 2f3i3 = 10 -5 , regardless of the viewing direction. Moreover, the event number 
for sin 2 1Q\^ = 10 -2 in the polar direction shows a much steeper shock effect than for 
sin 2 2f5i3 = 10~ 5 . Those features can be understood as follows. 

As already mentioned, the change in the event number is sensitive to the adiabaticity 



10" 



and there 



condition (e.g., Eq. (|2.4| ). If 613 is relatively large (e.g. sin 26*13 = 10~ 
is no shock wave in the resonance, the adiabatic condition is satisfied for the H-resonance, 
converting almost all v x to v e . Since the average energies of neutrinos in the supernova 
environment are in the following order, E Ue < E^ < E VxjDx , the average energy of v e 
becomes higher than the original z/ e due to the flavor conversion from v x . On the other 
hand, if #13 is small (e.g., 10 -5 ) or there is a shock wave in the resonance, the adiabatic 
condition is not satisfied. In this non-adiabatic case in the H-resonance, the original v e 
spectrum almost keeps unchanged. 



As stated in section 3.3, the cross section for detecting v e depends on E„. Therefore 



the v e event number of the adiabatic case becomes larger than that of the non-adiabatic 
case. If the resonance is closely non-adiabatic always, the shock effect becomes very minor, 



which is the case for the right panel in Figure 14. We also calculate R(H/L) denned by 
Eq.(|3,7j). In the case of sin 2 2#i3 = 1CP 2 , R(H/L) of u e is similar to Figure [l^. Conversely, 
the time evolution of R(H/L) of v e of sin 2 1Q\^ = 10 -5 does not show the shock feature. 

4.2 Difference of the original neutrino spectrum 

In this section, we discuss the variations of the original neutrino spectra, and its effect on 



the obtained results. As already mentioned in section 2.2, we have employed the one, which 
is assumed to decay exponentially with a timescale of ti, = 3 s with the average neutrino 
energies kept unchanged and is adjusted to be consistent with the total neutrino energies 
in the model of Lawrence Livermore group fi0|| . Such a treatment has been often employed 
so far in the supernova neutrino and nucleosynthesis studies (e.g. Yoshida et al. 2006 p3| 
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Figure 15: The evolution of average neu- 
trino energies for the time-dependent spectra 
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solid line (red line), dashed line (green line) 
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Figure 16: Same as Figure 15, but for the 
neutrino emission rate as a function of time 
inEq. (4.1) (e.g., 0). 



and reference therein). When we set ti, to be 6 sec, the obtained results are qualitatively 
unchanged though the event number becomes smaller up to a factor of ~ 2 than for the 
decay-time of 3 s. 

To explore the uncertainty of the original spectra, we here employ another original 
neutrino spectra used in [fiTf ] as follows, 



K = —jt<p( e » 



dt 



(4.1) 



where dN v /dt is a neutrino emission rate, and <p(E) is a normalized energy spectra defined 
as, 



<t>(E v 



r(a + i) {W)J (K) 



(4.2) 



(E), a, and T is the average energy, an energy shape parameter, and the gamma function, 



respectively |17]. For simplicity, we have taken a = 3 for all flavors [3C]. In this model 



^rp and (E) follow the calculation of the Lawrence Livermore group [1C]. The time evo- 



lution of the average energy and the neutrino emission rate is shown in Figure |15| and |16|, 
respectively. The solid line (red line), the dashed line (green line) and the dotted line (blue 
line) corresponding that of u e , v e and u x , respectively. Hereafter we call this spectra as 
"time-dependent" original spectra. 



Figure 17 shows the neutrino spectra at the SK calculated by the mentioned procedure 



in section 3.3, in which left and light panel is for v e and v e , respectively. It can be shown that 



the neutrino spectra are qualitatively insensitive to the difference in the original neutrino 
spectra examined here. However significant difference appears in the event number. Figure 



18 and 19 is the event number of v e and z^ e , respectively. The difference of the event 
number between the polar and the equatorial direction becomes much smaller than that 
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Left panels are spectra of D e , and right figures are v e . The solid lines (red lines) and the dotted 
lines (blue lines) are the spectra in the polar direction and the equatorial direction, respectively. 
We assume sin 2 2(9 13 = 10~ 3 and inverted mass hierarchy. 



of our fiducial model shown in Figures I and p]. This is because the neutrino flux of the 
time-dependent spectra decreases already before the shock reaches to the resonance layers. 

R(H/L) is shown in Figures [2(] (for v e ) and |2l] (for u e ). In the equatorial direction of 
Figure the ratio gradually increases although it is constant in Figure [!(]. This reflects 
the increase of the average energy of v e with time (see the dotted line of Figure |15| ). In 
contrast, the ratio of the equatorial direction is almost flat in Figure 21, reflecting the 
evolution of the average energy of v e . The error bars here are shown to be larger than 
for Figure [TJ[ simply because the event number of v e here becomes smaller than for the 
original neutrino spectrum in section 2.2 (compare Figure 12 with 19). Albeit sensitive to 
the feature of the original spectra, it is found that the ratio for the polar direction becomes 
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generally smaller than for the equatorial direction due to the shock effect (Figures 10, 20, 
21) for relatively large value of #13. 

Those results indicate that the shock effect on the event number discussed so far, is 
subject to the large uncertainty of the original neutrino spectra. What is qualitatively 
new pointed out in this paper is that in the MHD explosions, the timescale of the shock- 
arrival could be shorter than the decay timescale of the neutrino flux, albeit depending 
on the evolution of the neutrino spectra. It should be also mentioned that if the energy 
difference between v e and v x becomes small due to detailed neutrino interactions near the 
neutrinosphere [72, 73], the shock effect should be also small. To determine the original 
neutrino spectra in a consistent manner with the shock-dynamics in the supernova envelope 
requires more sophisticated MHD modeling of core-collapse supernova, in which detailed 
microphysics is also included. Several groups are really pursuing it with the use of advanced 
numerical techniques [46, 47], however still unfeasible at present. 

In this paper, we neglected the effect of neutrino self- interactions. Among a number of 
important effects possibly created by self- interactions (e.g., [29, |36|), we choose to consider 
the effect of single spectral swap of v e and v x suggested by Fogli et al. [74]. In our 
model, the original spectrum was cj)p e > <\> Vx in the low-energy side, and (pp e < (j> Vx in the 
high-energy side. Due to the spectral swap, the above inequality sign reserves. Thus the 
influence of the shock should occur also oppositely (see Eq.([0]) and discussion in section 
3.2). Therefore the event number of neutrinos in the polar direction is expected to increase 
by the influence of the shock. To draw a robust conclusion to the self-interaction effects 
on our models naturally requires more sophisticated analysis (e.g., |54j), which we pose as 
a next step of this study. 

As a prelude to more realistic computations of the flavor conversions taking into ac- 
count the effects of neutrino self-interaction, also with more realistic input of the neutrino 
spectra, our findings obtained in a very idealized modeling should be the very first step 
towards the prediction of the neutrino signals from magnetic explosions of massive stars. 



5. Summary 

We studied neutrino oscillations from core-collapse supernovae that produce magnetohy- 
drodynamic (MHD) explosions, which are attracting great attention recently as a possi- 
ble relevance to magnetars and/or long-duration gamma-ray bursts. Based on a recent 
supernova simulation producing MHD explosions till the shock break-out, we calculated 
numerically the flavor conversion in the highly non-spherical envelope through the pure 
matter-driven MSW effect. As for the neutrino spectra at the neutrino sphere which the 
MHD simulations lack, we employed the two variations based on the Lawrence Livermore 
simulation, with or without the evolution of the average neutrino energies. With these 
computations, we investigated how the explosion anisotropy could have impacts on the 
emergent neutrino spectra and also on the observed neutrino event number at the Super- 
Kamiokande detector for a Galactic supernova. 

In the case of the inverted mass hierarchy with a relatively large #13 (sin 2 26\^ > 10 -3 ), 
we demonstrated that survival probabilities of v e and v e seen along the rotational axis of 
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Figure 18: The expected event number of 
v e using the time-dependent original spec- 
tra. The solid lines (red lines) and the dotted 
lines (blue lines) are the event number for the 
polar and equatorial direction, respectively. 
The source is assumed to be at 10 kpc. 
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Figure 19: Same as Figure 18 but for v e 
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Figure 20: The ratio of event number of 
D e using the time-dependent original spectra. 
The solid lines (red lines) and the dotted lines 
(blue lines) are the event number in polar di- 
rection and equatorial direction, respectively. 
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Figure 21: Same as Figure 20 but for v e 



the MHD supernovae (i.e., polar direction), could be significantly different from those 
seen from the equatorial direction. The event numbers of v e observed from the polar 
direction show steepest decrease, reflecting the passage of the magneto-driven shock to the 
so-called high-resonance regions. We pointed out that such a shock effect, depending on the 
original neutrino spectra, appears also for the low-resonance regions, which could lead to a 
noticeable decrease in the v e signals. This reflects a unique nature of the magnetic explosion 
featuring a very early shock-arrival to the resonance regions, which is in sharp contrast to 
the delayed neutrino-driven supernova models in spherical symmetry. Our results suggested 
that the two features in the v e and v e signals, if visible to the Super-Kamiokande for a 
Galactic supernova, could be an observational signature of magneto-driven supernovae. 
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One of the nearest supernovae, associated with the long-duration gamma-ray burst 



(GRB) is, SN1998bw (Type Ic) at the distance of ~ 14 Mpc [58]. As mentioned earlier, 
such events are likely to be associated with the energetic stellar explosions induced by 
magnetic mechanisms, often referred to as hypernovae (e.g., J7j| and references therein). 
It is interesting to note that the fraction of SNe Ib/c associated with long GRBs is small 
(~ 1% — 10%), however, is appeared to coincide with that of hypernovae. Albeit with 
large uncertainties of GRB statistics, the link between hypernovae and long GRBs seems 
to be observationally supported (see discussion in and reference therein). Recently the 
proposals of Mton-class detectors such as Hyper-Kamiokande [[[(]] and Deep-TITAND {nj , 
are becoming a real possibility, by which the shock effect studied here could be visible out 
to the Megaparsec distance scales. This implies that the shock effect, if observed (much 
nicer if observed also by the electromagnetic and gravitational- wave observations |77], ^] ) , 
could provide important hints to understand the long-veiled central engines of gamma-ray 
bursts, because the shock arrival time to resonance layers should reflect the activities of 
the engine. We hope that this study could give strong momentum to theorists for making 
more precise predictions of the neutrino signals from magneto-driven supernovae. 
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